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1. Introduction 

Since the work of T.T. Wu et al. P|, it is known that the expansion in n-particle 
contributions to the zero field susceptibility of the square lattice Ising model at 
temperature T can be written as a sum: 



of {n — l)-dimensional integrals pi El IH HI EI 0) the sum being restricted to odd 
(respectively even) n for the high (respectively low) temperature case. While the first 
contribution in the sum, x'^\ is obtained directly without integration, and the second 
one, x^^\ is given in terms of elhptic integrals, no dosed forms for the higher order 
contributions are known despite the well-defined forms of these {n — l)-dimensional 
integrals. The x''"^'s are (n — l)-dimensional integrals of holonomic (algebraic) 
expressions, and are consequently holonomic, or "D- finite": they are solutions of finite 



oo 




(1) 
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order differential equations with polynomial coefficients [HI EI- Unfortunately, such 
theorems of holonomy theory do not give any lower or upper bound, or hint, as to the 
order of the differential equation satisfied by a given x'"-*. 

As far as singular points are concerned (physical or non-physical singularities in 
the complex plane), and besides the known s = ±1 and s = ±i singularities, B. Nickel 
showed jSl that 

^(2n+l) -g g-jj 

gular| for the following finite values of s = sh(2 J/kT) 
lying on the unit circle (m = k = excluded): 

2-(s + -) = + ^ + + — 

=1, -n < m, k < n (2) 

When n increases, the singularities of the higher-particle components of x(s) 
accumulate on the unit circle \s\ = 1. The existence of such a natural boundary 
for the total x(s), shows that x(s) is not D-finite (holonomic) as a function of s. 
To understand the analytical structure of such a transcendental function, it is thus 
crucial to better understand the analytical structure of the "holonomic" x'"^'s, and, 
as a first step, to find the still unknown differential equation verified by the 3-particle 
contribution x*''^''- 

A significant amount of work had already been performed to generate isotropic 
series coefficients for x^"^ (by B. Nickel 0d up to order 116, then to order 257 by A.J. 
Guttmann§ et al). More recently, W. Orrick et al. pHl, have generated coefficients|| of 
X(s) up to order 323 and 646 for high and low temperature series in s, using Perk's non- 
linear Painleve difference equations for the correlation functions |1()[ IllL fT^ I13L I14j . 
As a consequence of this non-linear Painleve difference equation and the associated 
remarkable recursion on the coefficients, the computer algorithm had a 0{N^) 
polynomial growth of the calculation of the series expansion instead of an exponential 
growth that one would expect at first sight. However, in such a non-Hnear, non- 
holonomic, Painleve-oriented approach, one obtains results directly for the total 
susceptibility x{s) which does not satisfy any linear differential equation, and thus 
prevents the easily disentangling of the contributions of the various holonomic x^"'''s. 

In contrast, we develop here, a strictly holonomic approach. This approach 
enabled us to get 490 coefficients of the series expansion of x'"^\ from which we 
have deduced the Fuchsian differential equation of order seven satisfied by x^'^^- 
This Fuchsian differential equation presents a large set of remarkable properties and 
structures that are briefiy sketched here and will be analyzed in more details in 
forthcoming pubHcations, along with the analytical behavior of the solutions. The 
method used in this paper to obtain the Fuchsian differential equation is not specific 
of the third contribution x'''^'' and can be generalized, mutatis mutandis, to the other 
x'"' , (n > 3) without any drastic changes in the mathematical framework^. 

I The singularities being logarithmic branch points of order ^2n(n+i)-i . with e = 1 — s/si 
where is one of the solutions of |2j. 

§ Private communication. 

II The short-distance terms were shown to have the form (T — Tc)'' ■ {log\T — Tc|)' with p> ■ 
% But certainly requiring larger computer calculations. 
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2. Generating the series for x^^^ 

Similarly to Nickel's papers [£,^2,, we start using the multiple integral form of the 
X^"^'s. We focus on x'^^ and consider the double integral: 

(^ «4\i/4 



S 

(.) = ^ r r '^^^ y- h^^' m 

47r^ 7o Jo \l-XiX2X'iJ 
Xj = ^ , (4) 

1 + — S cos (pj + ^y {1 + s'^ — S COS — 

y^ = rrr-^ i ^ x ^ l=^' j = 1,2,3, 0i + 02 + 03 = 

\J\ \ + — scosqyjY — 

Many forms [B, 7^ for iJ^^^ may be taken and are equivalent for integration purposes, 
eg-, 

H^'^ = /23 (/ai + = (sin0. - sin0,) (5) 

It is straightforward to see that x'''^'' is only a function of the variable w = ^s/{l + s'^). 
From now on, we thus focus on x^'^^ seen as a function of the well-suited variable w 
instead of s. We may expand the integrand in ijSl in this variable w and integrate 
the angular part. For x'''^'', this would mean there are 18 sums to carry out. We, 
instead, expand the integrand in the variable x oi where we succeeded in deriving 
remarkable formulas for yi" carrying one summation index. As a consequence, 
we are able to write x^^^w) as a fully integrated expansion. Our algorithm runs 
in a polynomial time calculations (namely 0{N'^)). In contrast with Orrick et al. 
calculations ^H], this calculation is not based on any recursion: it allows one to 
obtain any given coefHcient separately without requiring the storage of all the previous 
data. The details of these calculations, tricks+, and of this program will be given 
elsewhere |T5] . 

At present, we have obtained* the expansion of x''^'' up to w'^^^ . As expected, 
this expansion is in agreement with the previous results published by Nickel j^l C] 
using a numerical method of integration, as well as the improved unpublished results 
(up to w^^"^) by A.J. Guttmann U et al. In terms of the well-suited variable w, the 
first terms of the expansion of x'^^ {w) read : 

^i^^ = w'3-H36w"-H4wi2^_884y^i3^ j^9g^i4 
8 

+ 18532 w^^ 6084 w^^ + ■ ■ ■ (6) 



with: 



3. The Fuchsian differential equation satisfied by x^^\w) 

Given the expansion of x'^H""^) up to the next step will be to encode all the 

numbers in this long series into a linear differential equation. Note that such an 
equation should exist iSii^I though, its order is unknown. Using a dedicated program 

+ For instance, our calculations also underline the important role played by hypergeometric functions. 
* We can get a new coefficient every two days, 
(t Private communication. 
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for searching for such a hnear differential equation with polynomial coefficients in w 
and steadily increasing the order, we succeeded finally in finding the following Hnear 
differential equation of order seven satisfied by the 490 terms we have calculated for 



with : 



n— 



:F{w) = (7) 



a? = 

w"^ -{l-w) (l + 2w) {l-Awf{l + Awf (l + 3w + 4w2) Pj{w) 
ae = w'^ ■ {I ~ Awf {1 + 4,wf Pe{w) 
az = ■ (l-4w)^ (l + 4w) Pr^{w) 

a4 = . (1 -Awf Pi{w), ag = -{l-Aw) P^iw) 

02 = V? P2{w), ai = wPi(w), flo = Pn{w) (8) 

where Pt{w),Pq{'w) ■ ■ ■, Po{w) are polynomials of degree respectively 28, 34, 36, 38, 
39, 40, 40 and 36 in w p^I- These polynomials are too large to be given here. 

Note that the series of x'"^\w) up to order 359, is sufficient to identify 
the differential equation when the order q = 7 and the successive degrees 
47, 46, 45, 44, 43, 42, 41 and 36 of the polynomials in front of the F(")(w)'s are imposed. 
Inside this framework, and since our series has 490 coefficients, we have here 131 
verifications of the correctness of this differential equation. 

At this point, let us remark that if one had worked with variable s, instead of w, 
the series in s up to order 699 would be needed in order to obtain the seventh order 
differential equation satisfied by x'^'^H'*)- 

Since the singular points of this differential equation correspond to the roots of 
the polynomial corresponding to the highest order derivative namely 07, we 

give the exact expression of P7{w) : 

P^(w) = 1568+ 1563810- 565286 - 276893^3 
+ 34839063 + 100696470 - 1203580072 

- 5514282112 18005067728 110343422816 

- 140604884224^1°- 1825536178688 u;" + 920432273408 
+ 28913052344320 u;" + 38181758402560 w^'' 

- 112307319603200 w^^- 544140071665664 u;^^ 

- 1144172108054528^1^ - 1027222993371136^1^ 

- 1992177026596864 - 2948885085421568 
+ 221 1524294737920 + 8204389336481792 

+ 675795924156416 - 2882636020187136 «724 

- 5364860829302784 - 222238787764224^2'^ 

+ 158329674399744 + 39582418599936 (9) 

The differential equation Q is an equation of the Fuchsian type since there are 
no singular points, finite or infinite, other than regular singular points. With this 
property, using Frobenius method ^^I; it is straightforward to obtain, from the indicial 
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equation, the critical exponents, in w, for each regular singular point. These are given 
in Table 1. 





s = 






9,3,2,2,1,1,1 


w = -1/4 


s = -1 




P = 


3,2,1,0,0,0,-1/2 


w = 1/4 


s = 1 




P = 


1,0,0,0,-1,-1,-3/2 


w = -1/2 


1 + S + s2 ^ 




P = 


5,4,3,3,2,1,0 


w — 1 


2 - s + 2s2 = 




P = 


5,4,3,3,2,1,0 


l + 3w + 4w2 = 


(2s2 + s + l)(s2 ^ 5 _ 


h2) = 


P = 


5,4,3,2,1,1,0 


l/w = 


1 + s2 = 




P = 


3,2,1,1,1,0,0 


w = wp, 28 roots 


s — sp, 56 roots 




P = 


7,5,4,3,2,1,0 



(m- 1) g (g- 1) 



(10) 



Table 1: Critical exponents for each regular singular point, wp is any of the 28 
roots of P7{w). We have also shown the corresponding roots in the s variableff. 

At this point, it is worth recalling the Fuchsian relation on Fuchsian type 
equations. Denoting wi, W2, ■ ■ ■ , Wm, Wm+i = oo, the regular singular points 
of a Fuchsian type equation of order q and pj.i, • • • , pj^q ( j = 1, • • • , m + 1) the 
q roots of the indicial equation corresponding to each regular singular point Wj , the 
following Fuchsian relation holds : 

m+l q 

H H ^J'fe 
i=l k=l 

The number of regular singular points is here m+l = 36 corresponding respectively 
to the 28 roots of P7, the two roots of 1 + 310 + iw^, five regular singular points, 
and the point at infinity w — 00. The Fuchsian relation is actually verified here with 
q = 7, m — 35 , the sum of all the Pj,kS being actually 714 . 

Let us comment on the singularities appearing in 10. The well-known 
ferromagnetic and antiferromagnetic critical points correspond respectively to w = 
1/4 and w = —1/4, while w = corresponds to the zero or infinite temperature 
points. The value w = 00, or s ~ ±i, is known as a non-physical singularity. 
The values w = 1 and w = —1/2 correspond to the non-physical singularities ltTT| 
described by Nickel for n — 1 {m = k — excluded) : 



1 k 
— = 

w 



1 



^ 1, 



1 



-1 < m, fc < 1 



(11) 



F\irthermore, besides the known singularities mentioned above, we remark the 
occurrence of the roots of the polynomial P7 of degree 28 in and the two quadratic 

ft In the variable s the local exponents for w = ±1/4 are twice those given. 
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numbers l + 3u; + 4u;^ — which are not of the form ((TTll . The two quadratic 
numbers are not on the s-unit circle : \s\ = \/2 and |s| = 1/a/2. 

To analyze the local solutions of the differential equation, let us recall that, in 
general, it is known QS] that for a set of k local exponents such that the difference, in 
absolute value, between any two of them is an integer, the local solutions may contain 
logarithmic terms up to log*^"^. In fact, the Fuchsian equation ijjjl is such that the 
occurrence of logarithmic terms, in the local solutions near a given regular singular 
point, is due only to the occurrence of multiple roots in the corresponding indicial 
equation: a root of multiplicity p inducing logarithmic terms up to log''"^, i.e., at 
most log^ for Q , as it is shown also in the monodromy matrices (see below) . 

More precisely for P-j, near any of its roots, all the local solutions carry no 
logarithmic terms and are analytical since the exponents are all positive integers. The 
roots of P7 are thus apparent singularities of the Fuchsian equation Q. Note 
that the "apparent" character of P7, means that, to ensure the absence of logarithmic 
terms in the general solution, there are (g — l)/2 = 21 relations between the various 
P„'s at each root of P-j (for more details see flE!). 

The two unexpected quadratic numbers correspond to singularities of solutions 
of which behave locally like (1 + 3w +4 vP') ■ ln(l + + 4 w^) ■ qi +92, where qi 
and 92 are analytic functions near the two quadratic roots 1 + 3w + 4 = 0. This 
"weakly" singular behavior can also be seen on the monodromy matrix (see below) 
associated with these two roots, where one finds a nilpotent matrix of order two (no 
log^ term). 

Details on the local solutions of around each regular singular point, together 
with the constants giving the particular physical solution x^'^\ will be given in a 
forthcoming publication. 

We will sketch here the analysis of the Fuchsian equation by focusing on the 
local solutions of Q around w = 0. Recall that the exponents near w = Q are 
respectively: 9, 3, 2,1 and again 2,1, 1, all of them integers. Besides the x'-'^V^ 
solution of the Fuchsian differential equation which behaves, near w = Hke lO, 
(which we will denote, from now on, Sg, since its leading term behaves like w^), one 
expects solutions with leading terms behaving like w , and w'^ , because of the 
previous integer exponents near w — 0. Furthermore, because of the singularity 
confluence (repetition of same exponents) and the occurrence of integer exponents, 
one also expects solutions with the above-mentioned logarithmic terms. 

Actually, we have found two remarkable rational and algebraic solutions of ((Jjl, 
namely : 

^1 = S, = (12) 

l-4w' (l-4u>)Vl- 16w2 ^ ' 

a solution behaving like w'^ , that we denote S'3 : 

5*3 = w^ + 3w'' + 22w^ + 74w^ + 417w^ + 1465w** 

+ 7479 + 26839 w^^ + ■ ■ ■ (13) 

(2) (2) (31 

and three solutions S{ , S2 and S{ with logarithmic terms, behaving, at small w, 
as: 



S\^> = w ■ ln(w) • ci - 32 10* • ca (14) 
S'i'^ = ln(w) • di + 8 • da (15) 
Sf^ = 3 w • \n{wf • ei - 12 • ln(w) -62 - 19 ix;^ • 63 (16) 
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where ci, C2, di, ^2, ei, 62, 63 denote functions that are analytical at w = 0. 
However, these functions are not all independent. They can also be written in terms 
of S\, Si, S3 and as follows : 



SI 

s. 



(2) 



ln(w) • (S*! - 4 S'2 + 16 53 - 216 Sg) - 32 • C2 
ln{w) -{Si -2 S3 + 24 59) + 8w^-d2 



(17) 



(2) _ 
2 ~ 

5*^ ^ = 3 ln(w)2 . {Si +552-253) 

:.(2) 



6 ln(w;) • (2 53 - S\"' - 9 5^ ^) -19w''- 63 



where 



C2 

63 



w - 



167 

"96" 

5 

- w4 
7693 



2273 



6977 



96 



■ w 



120 



103 



315 



575593 



+ 

5 2191 

+ — 
2561473 



19371 



■ w 



40 

,4 I 



127434803 



456 11400 5700 93100 

Denoting Q — 2 in, one immediately deduces from the previous relations ifTTjl . the 
monodromy around w = : 



SI 



(2) 



SI 



(2) 



5 



(2) 



5 



(2) 



- n-{Si 

n-{S2' 



-452 + 1653 
253 +2459) 



216 5q 



5[ 



(3) 



5[ 



(3) 



(2) 



9 5 



(2)^ 



6n-{2S3 -S[ 

hsn^-iSi +552 -253) 



This calculation is a straight consequence of the fact that one only has a "logarithmic" 
monodromy, which just amounts to changing ln(w) into ln(w) + fl, in the previous 
expressions. Denoting Idr the 7x7 identity matrix, and taking the following order 
5i, 52, 53, 59, 5^^\ S^\ 5p'' for our seven-dimensional basis of solutions of 
one can see that any n-th power of the 7x7 monodromy matrix M, M„ = Af", 
satisfies — 3M^ + 3 Af„ — Idr = 0, and can be written as the following sum of 
the identity matrix, and the two first powers of the order-three nilpotent matrix N 
{N^ = 0) : 



= idj + 71 ■ n ■ N 


(n ■ n 

+ 2 


iV)2 




exp(n 


n-N) 








" 











1 
























-4 


1 





















16 


-2 - 


-12 


where : 
















-216 


24 



























6 
























54 




























(18) 



As it should, one can check that the form ifTsll . valid for any positive or negative 
integer n, is such that : M" • Af™ = 
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The monodromy matrices corresponding to the other regular singular points yield 
very similar results. Those corresponding to w — ±1/4 are matrices of determinant 
— 1, while all the others are of determinant +1. The monodromy around the 28 
roots of P-j is of course trivial (apparent singularities). The monodromy matrices 
corresponding to the other regular singular points verify very simple relations like the 
previous relation Ml - iMl + 3M„ - Idj = or Ml - 2M„ + /dy = 0, and 
also Ml - SMl + 3Ml - Idj = for those corresponding to w = ±1/4. This 
is a straight consequence of the fact that, in a Jordan form, only very simple Jordan 
blocks occurf, such as : 



J2 = 



(19) 



For instance, the monodromy matrices around the points w = and w ^ co are on 
the same footing, both requiring one J3 Jordan block, one J2 block, and the 2x2 
identity matrix. There is another set of monodromy matrices corresponding to the 
points w = 1, w = —1/2 and the two roots ofl+3w+4w^ = (associated 
with J2 blocks). More precisely, the monodromy matrices for w = and w = co 
verify Ml - 3 Ml + 3M„ - /dy = 0, but for w = -1/2 and w = 1 they verify 
Ml -2Mn + Idj = 0, as well as for the two roots of 1 + 3 w + 4 = 0. 

A detailed analysis of the monodromy group will be performed elsewhere. 
However, such an analysis of the Galois group| becomes easier to perform when taking 
into account some remarkable factorization and decomposition properties that are 
sketched in the next section. 



4. Algebraic properties of the Fuchsian equation 

Let us define : 

A = 

{l + 3w + Aw^f (l + 2u.)^ (-l + 4w)^^/' + (l - wf 

/ = (1 - w) (1 + 2w) (1 - Awf (1 + 4w)^ (1 + 3w + 4w2) 

We denote L7 the seventh order linear differential operator, corresponding to the 
Fuchsian differential equation 

where the QkiLiYs are the polynomials at defined in JHJ. Polynomial a'j{L'j) reads 
a7{L-!) = w"^ ■ f ■ P7, where we distinguish the "actual" and the apparent singularities. 

We now give some remarkable factorization properties of the linear differential 
operator L7. The fact that the very simple expression Si (see Ijl2ll l is a solution of 
the differential operator Lj and that it is also solution of the first order differential 
operator Li : 



duu w {1 — iw ) 
t We thank Jacques-Arthur Weil for this result. 

I The main difficulty is to find a global structure like the monodromy Galois group from the knowledge 
of all these local monodromy matrices expressed in the different well-suited local basis associated with 
each regular singular point. 
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implies the following factorization of L7 (or more precisely the right-division of L7 by 
Lx): 

Lr = Me ■ Li (22) 

Similarly, the fact that S2 is a solution of ^7, and that it is also clearly solution of a 
first order differential operator A^i : 

^ A _ 2(l + 2z^) 
dw (1 — 16 w^) 
implies the existence of another factorization of L7 : 

L7 = Nfi- Ni (24) 
The linear differential operators of order 6 read {Xe denoting Mg or iVe) : 



^6 1 5 J; 

where ak{Xe) are polynomials in w and a6(Af6) = a6(^6) = ""^^ • / • Priw). 

These two factorizations are consequences of the existence of remarkable simple 
algebraic solutions of L7. Besides this, there exists another factorization of L7 related 
to the adjoint differential equation of (3.1). This is explained as follows. One can also 
see that the adjoint of L7, denoted L7, admits the following rational solution§: 

SI ^ (26) 

with : 

Qe = 1 + 19 w - 368 - 3296 + 17882 + 272599 

+ 160900 - 6979208 + 7550800 + 203094872 

- 278920192 - 3959814304 - 2115447424 w^^ 

+ 20894729472^13 + 39719728128 w^'' + 20516098048 

+ 256763363328 w^^ - 327065010176 - 8810227761152 w^^ 

+ 414933057536 w^^ + 116411936538624 w^o 

+ 296827723186176 w^^ + 317648030138368 w^^ 

+ 179148186189824 + 194933533179904 u;^^ 

+ 112931870081024^2^ - 55246164328448 w^'^ 

+ 11063835754496 + 1511828488192 w^^ (27) 

yielding immediately the following factorization (or more precisely the left- 
division of Li by Ml): 

L7 = Mi-Lg, where: Lq = L5 ■ Ni (28) 

d 1 dSl , , 

with: Ml = -r+ sz-r- 29 

dw S\ dw 

and where L5 and are linear differential operators of order five and six respectively, 
which read : 

di 1 '"""^ (t 

§ We thank Jacques-Arthur Weil for the remarkable result 1261 and 1281 . 
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with : 

06(^6) = ■ f ■ Qeiw) 

a^iLe) = w"^ ■ {1 - 4w)^ {1 + Awf ■ Q5H (31) 

05(^5) = (l-4u;)(l + 4u;)'-/-Q6H 

The roots of polynomial Qe (w) axe apparent singularities of the differential equations 
associated to the operators Lq and L5. 

One can also see that Mg in Ij22ll can also be decomposed as follows : 

Me = Ms • Ti with: 



where a^iM^) ^ {1 - Awf (1 + Awf ■ f ■ Priw), and: 

^ _ d l + 4w + 48w^ 
^ dw w (1 — 16 w^) 

Calculations seem to show that L5 (resp. M5) is irreducible : we have not found any 
further factorizations like L5 — L4 ■ A\ or L5 ~ Bi ■ L4, or even L5 = £3 • A2 or 
L5 — B2 ■ L3. 

One can also be interested in the Wronskians of L7 and Lq. They can be written 
as follows : 

WiLr) = ^, WiLe) . (32) 

thus checking the following relation on the Wronskians of L7, Lq, and Mi, deduced 
from Ij28ll . namely: 

W^(L7) = W^(Mi) • VF(i6) (33) 

From these factorizations, and/or decompositions, relations, corresponding to the 
existence of rational, or algebraic (square root of rational), expressions, one should not 
be surprised to find that the Wronskians of ah the various operators we define here, 
are remarkable rational, or algebraic (square root of rational), expressions (see l(32)l ). 
The fact that the squares of the Wronskians ll32ll are rational functions corresponds 
to the following identities between P7 and Pg (resp. Qe and Q^) : 

A' (l-4w)4(l + 4w)2 Pe 



-jk 



P7 A wf Pt 



= (34) 



g^-i+4-4+"-'"°''''-^^°'''.gi = ,35) 

Q& w f \ wf Qe 

The Wronskians (or their inverse for the adjoints) of ah the differential operators 
L7, L7, M6,Mg, M5, M5 can be expressed as simple powers of P7 divided by a 
simple expression similar to A in lf2?Hl . Similarly, the Wronskians (or their inverse 
for the adjoints) of all the differential operators Lq, Lq, L5, Lg can be expressed as 
simple powers of Qq divided by a simple A-like expression. All these Wronskians are 
thus rational (L5, Lg, M5, M5) or such that their square is rational (L7, L7, Lq, 
Lq, Mq, Mg). Of course, there are relations between these Wronskians which are in 
agreement with the operator factorizations previously described (see l|33jl V 

Besides the "fundamental" singularities w = 0, —1/2, ±1/4, 1 and the 
"unexpected" quadratic numbers singularities solutions of 1 + + A-uP' = 0, one 
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could say that a large part of the "arithmetic complexity" of the operator L-j arises 
from the two polynomials of degree 28 in w, Pr and Qq. 

The occurrence of the apparent singularities, associated with the quite large 
polynomial P7, can be considered slightly unpleasant or disturbing : one would like 
to exchange the seventh order Fuchsian equation , for another differential equation 
(or a differential system) where these "spurious" singularities have disappeared. This 
is the so-called desingularization problem for differential equations ^HilSl- We have 
performed such "desingularization" : the "price" to be paid is that one has no longer 
a seventh order differential equation, but an eighth order differential equation. The 
associated eighth order differential operator reads: 

1 dJ' 

{l-AwY^ (l + 4w)"^ {\ + 2wf' (l-w)"= (l + 3w + 4u;2)"« 

ni, • • • , ng being positive integers. The other polynomials ak{L%), are quite large and 
will not be given here. Let us sketch this desingularization procedure. Recall the 
exponents associated with the roots of P-j (see Table 1). They are 0,1,2,3,4,5,7; 
one sees that the exponent 6 is missing. Basically, the method amounts to building 
a differential equation having all the solutions of 10, together with a solutionff 
associated with this missing exponent 6, in order to "fill the gap" and make the roots 
of P-j ordinary points for the homogenous differential equation associated with L^. 
Note too that the desingularized differential equation, or equivalently the operator Lg 
are far from being unique. The various eighth order Lg operators are, by construction, 
of the form Lg = Li • L7, where the first order operator Li is quite involved^. Prom 
a desingularized form Hke l(3fijl one can also introduce a differential system+ : 

e-Y ^ A-Y, Y = (y, 0-y, 0^ -y, 0^ -y, ■■■) (37) 

-,1, , Ao , Ai , A2 .A3 

with : A = — 



w w-l/4: w + 1/4 u; + l/2 
+ + — T. ; ;— + P(w) ■ Ids 

where /dg denotes the 8x8 identity matrix, P{w) denotes a polynomial in w, and 
where the A„ matrices are very simple 8x8 matrices, simply related to the monodromy 
matrices, and where denotes the "well-suited" derivation associated with the "true" 
singularities : 

^ w{l-w) {l + 2w) (l-4u;) (l-|-4w) il + iw + Av/) ~ 

aw 

A form like Ij37|l is clearly much simpler, and more canonical : the calculations one has 
in mind (Galois group, rigidity index (see below)) should be much simpler to perform 
with this canonical system form l|37ll . Details will be given elsewhere. 

ft The solution 5 = is too naive, it puicks out the apparent singularities of P7 but introduces 
new spurious apparent singularities. 

1 As far as the solutions of Lg are concerned, the differential operator Lg adds to the known solutions 
of -L7, a solution with a strong exponential behavior. 
+ See for instance the chapter 6 of 
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5. Comments and Speculations 

A better understanding of the total susceptibility x certainly requires an exhaustive 
knowledge of the singularities of the successive Fuchsian differential equations 
associated with the n-particle contributions x^"'*- Besides the apparent singularities 
associated with the roots of the polynomial P7, we have noted the occurrence of the 
two rather unexpected quadratic numbers solutions of 1 + 3 w + 4 w'^ — 0, which, also, 
are not of Nickel's form lfTT|l . The eUiptic parameterization of the Onsager model is 
well-known. Recalling the exact expression of the modular invariant M for the Ising 
model (see for instance (201) : 



^=7^ ....,2/ (38) 



1 (1- 16^2 + 16^4 

1728 w8 (1 - 16^2) 

it can be seen that the two quadratic numbers correspond to a rational value of the 
modular invariant: Ai = —125/64, while w = 1 also corresponds to a rational 
value M — —1/25920, that w — ±1/4, 0, 00 correspond to A4 = 00, and that 
w = —1/2 corresponds to the rational value A4 = 32/81. Could this mean that 
the Fuchsian equation 10 could be "canonically" associated with an elliptic curve ? 
One can, for instance, recall the period mappings and Picard-Fuchs equations (and 
beyond mirror symmetries ...) associated with a family of elliptic curves |21II22) . The 
Picard-Fuchs equation : 

144s (s- 1)2. (s-^ + ^)+ (31* -4)-/ = (39) 
is associated with the family of elliptic curves : 

2 . s 27 27s 

= 4 x"^ H • X H 

1 — s 1 — s 

We now make a few comments on the physical solution Sg. The Fuchsian 
equation is highly non-trivial and structured : from the previous factorizations, 
and/or decompositions, one might imagine that 5*9 — x^^V^ is in fact a solution 
of a sixth order homogeneous differential equation, or even a fifth order homogeneous 
differential equation. This is not the case. Relation ll28|l means that 5*9 — x^'^V^ 
can be seen as a linear combination of Si and of a solution of a sixth order linear 
homogeneous differential equation, S{Lq), associated with the operator Lg. The 
solution 5*9 = x^^V^ is actually such a linear combination Sg = a-Si+ S{Lq), a 
being different from zero*. The three-particle contribution, x'''^'', is thus a solution of 
the seventh order differential equation Q and not of a homogeneous linear differential 
of smaller order. 

Coming back to the analysis of the Galois monodromy group [231 I24L 
of the Fuchsian equation it is clear that the previous factorizations, and/or 
decompositions, impose severe reductions of the group. From these various 
factorizations, or decompositions, of the linear operator Lj = Mi-Lg — Mi- L^- Ni, 
one can deduce that the Galois monodromy group is isomorphic to Gal{LQ), the Galois 
monodromy group of Lg- Introducing Gal{L^), the Galois monodromy group of L^, 
one deduces from Lq = L5 • Ni, that the operator Ni injects Gal(Lc^) in Gal^Lo). 
This does not mean that Gal{Le) is isomorphic to Gal{Lz), but that knowledge of 
Gal[L^) is required to describe Gal{Le). RecaUing the rationality of the Wronskian 

* The coefRcient a characterizing the "projection" of Sg on Si, can easily be calculated writing 
LeiSg - a- Si) = 0. One finds a = 1/24. 
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of L5, one deduces that Gal{Lr,) is a subgroup of SL{5,C) and not GL(5, C) (the 
rationality of the Wronskian means that all the monodromy matrices of L5 have +1 
determinants). The monodromy matrices of Lq (as well as the one of the Fuchsian 
equation Q) have ±1 determinants. The Galois group Gal{Le) is, up to a Z2- 
graduation, a subgroup of SL{6,C). Therefore the Galois group of the Fuchsian 
equation ijjjl is represented by 6 x 6 matrices with ±1 determinants, most of its 
structure requiring the analysis of a Gal{L^) subgroup of SL{b,C). Of course, one 
can also understand the Galois group of the Fuchsian equation 10 from the analysis 
of the Galois group of M5, namely Gal{M^) which is also subgroup of SL{5,G) 
(isomorphic to Gal{L^)). It seems that L5 is not self-adjointft (modulo conjugation 
by an operator) and irreducible (see Section 4). The two facts seem to rule out the 
existence of a symplectic structure. A more detailed analysis of the Galois group of 
the Fuchsian equation 10, that is, of the Galois groups Gal{L5), and Gal{Le), will 
be given in forthcoming publications. 

Let us, now, focus again on the "physical" solution of the Fuchsian equation 
'S'g = x'"^V8) and on its successive derivatives Sg, Sg , Sq^\ - ■ ■ We have seen that, 
as far as the linear dependence of these expressions is concerned, we have a seven- 
dimensional vector space {Sg is not a solution of a homogeneous sixth order differential 
equationf). However, to better understand the "true nature" of the susceptibility Xj 
one would like to characterize the "degree of transcendence" of Sg = x^^V^) that 
is the minimal number of successive derivatives of Sg satisfying an algebraic non- 
linear relation"'". The Galois monodromy group gives a valuable information about 
this "degree of transcendence". Let us consider the orbit of 6*9 under the Galois 
monodromy group: one gets an algebraic variety dense in the subspace of solutions. 
The dimension of this algebraic variety is actually this degree of transcendence. This 
analysis, however, requires an exact knowledge of the Galois monodromy group. 

The susceptibility x has been shown to be a transcendental (non-holonomic, non 
D-finite) function: it cannot be solution of a linear differential equation, but this does 
not mean that it cannot be solution of a differential equation. Along the previous 
"non-linear" line, one should emphasize that the possibility that x could be solution 
of some (Painleve-like ?) non-linear differential equation, is not yet totally ruled out ! 

Within the "linear" Picard-Vessiot Galois monodromy group framework, one can 
also try to evaluate the "index of rigidity j22l|23" of our Fuchsian differential equation, 
or the "index of rigidity" of the operators Lq or L5. Roughly speaking, this index of 
rigidity corresponds to the number of parameters that can be introduced to deform 
the linear differential equation, keeping the local monodromy matrices fixed|. For 
instance, the hypergeometric functions are totally rigid : as a consequence, "almost 
everything" can be calculated on such functions§. The differential operators L7, Lg or 
I/5 have many integer exponents, and, thus, their solutions have logarithmic behaviors 
{w — Wo)" • {ln{w — wq)) around each regular singular point wq. However, recalling 

tt However it can be shown that Le is not self-adjoint (modulo conjugation by an operator). La has 
S2 as a solution. If it were self-adjoint, Lg would have a similar quadratic solution which is not the 
case. 

t To be totally rigorous one should add that the minimal operator of Sg is a factor of L7, and since 
these factors are of order 1 or 6, and that Sg does not vanish on these factors, it must vanish on -L7. 
+ For instance the Weierstrass V function verifies the non-linear relation V'^ = iV'^ + g2'P + gs. 
X See the notion of rigid local systems 1271 1281 . Basically it amounts to calculating the dimensions of 
the centralizers of all the monodromy matrices. 

§ In the case of totally rigid systems, N. Katz has shown that the solutions have a geometric 
interpretation : they can be seen as periods of some algebraic varieties |27| . 
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idJ, itTfill and ltT6|l . and, more generally, considering the behavior near all the singular 
points, one sees that only Iniw — wq) and (ln{w — wq)) behaviors occur*, thus 
downsizing (see for instance (4.31) in Chapter (4) of j^) the index of rigidity of 
these differential equations. This is a strong indication that the Fuchsian differential 
equation Q, or the differential equations associated with Lq or L5, are extremely 
rigid. This can also be seen on the various monodromy matrices. The calculation 
of the index of rigidity, or in other words, the calculation of the "small" number of 
deformation parameters of ^ (resp. Lq, L^) will be performed elsewhere. This 
remarkable rigidity is not a surprise, when recalling the well-known isomonodromy 
theory of the Ising model, and in particular the occurrence of Painleve equations 
for the correlation functions of the Ising model. Does the deformation theory of 
10 in this "small" number of deformation parameters actually yield Painleve-like 
equations ? This is an open question. From a more down-to-earth viewpoint, this 
rigidity can be seen as "inherited" from the "total" rigidity of hypergeometric functions: 
our calculations for generating the w-series of x'^^^ are actually "flooded" by Hnear 
combinations (with binomial coefficients) of products of hypergeometric functions (see 
1 15) for more details). 

6. Conclusion 

The Hnear differential equation we have found for the three-particle contribution, 
X^^-* (''*)) to the susceptibility of the square lattice Ising model, is a highly structured 
and remarkable Fuschsian differential equation. We have sketched many of its 
remarkable properties and symmetries. It is also worth recalling that the 28 roots 
of the polynomial P7 are apparent singularities. We have been able to desingularise 
the Fuchsian equation ijjjl , in order to get rid of these apparent singularities related 
to the quite involved polynomial Pj. A deeper analysis of this Fuchsian equation 
(monodromy group, critical behavior of the solutions around the various singular 
points, ... ) will be given in forthcoming publications. 

All these results can be generalized, mutatis mutandis, to deduce the Fuchsian 
equations corresponding to the other n-particle contributions x'"'"s- The building of 
a computer program with a polynomial growth algorithm which can be generalized, 
mutatis mutandis, for the other x*^"^'s, was the key ingredient to get our Fuchsian 
equation. The ideas developed to create such polynomial growth programs underline 
the role played by hypergeometric functions, coming from a large number of remarkable 
identities on the underlying variables of the problem. One may also think that quite 
compHcated "fusion- type" relations on these hypergeometric functions can exist. This 
is crystal clear in the case of the 2-particle contribution x*"^^ dl- This however 
remains to be done in full. 

Clearly, beyond x^^'' ; a global understanding of the structure of the hierarchy of 
all the x'"^ 's could be contemplated. A better understanding of the total susceptibility 
X certainly requires an exhaustive knowledge of the singularities of the successive n- 
particle contributions x'"'' j or equivalently, of the corresponding successive Fuchsian 
equations. As far as analytical properties are concerned, we saw the occurrence in 
10 of unexpected quadratic numbers singularities, 1 -I- 3w -I- Aw'^ — (these two 
singularities are not on the s-unit circle: \s\ — \/2, Curiously these unexpected 

singularities correspond to a rational value of the modular invariant. Could this mean 

* This can be seen directly on the very simple Jordan blocks for the monodromy matrices (see 1131 '!. 
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that the Fuchsian equations for the successive x 's are "canonically^" associated with 
an eUiptic curve ? Should we rather understand this hierarchy of Fuchsian equations 
in a Hnear (or non-Unear) monodromy deformation theory framework ? 
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